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Sodium cobaltate at high Na content is a promising thermoelectric material. It has been recently 
shown that oxygen vacancies can alter the material properties, reducing its figure of merit. However, 
experimental data about the oxygen stoichiometry are contradictory. We therefore studied the 
formation of oxygen vacancies in Na^OoC^ with first principles calculations, focusing on x=0.75. 
We show that a very low oxygen vacancy concentration is expected at the temperatures and partial 
pressures relevant for applications. 



In recent years layered sodium cobaltate, Na x Co02, 
has attracted much attention because of its remarkable 
electronic and magnetic properties. Most of the interest- 
ing physics of this material comes from its quasi-2 dimen- 
sional cobalt oxide layers, surrounded by intercalating Na 
ions[l-3j. By changing the sodium content, one can con- 
tinuously tune the oxidation state of the cobalt ions. At 
high Na content, the majority of Co ions are in a 3+ 
oxidation state, and charge carriers are diluted holes. In 
this regime, the thermopower is particularly high prob- 
ably due to the hole hopping mechanism occurring from 
Co 4+ to Co 3+ along the oxide planes [4j [5] . 
The influence of point defects other than Na vacancies 
on sodium cobaltate is still debated. It has been re- 
cently reported that oxygen vacancies could have an in- 
teresting two-fold effect on the thermoelectric properties. 
Firstly, they can affect the cobalt oxidation state, i.e. 
reducing charge carrier concentration and increasing the 
thermopower. Unfortunately, such defects also lower the 
phonon scattering so that the overall effect is a reduction 
of the figure of merit of Na^CoC^-^ with increasing the 
defects concentration (<$) |6J. 

Experimental data about the concentration of vacan- 
cies in sodium cobaltate are contradictory, especially at 
x ~ 0.75. While some groups have identified an oxy- 
gen non-stoichiometry up to S = 0.16 [6-9J, others found 
no vacancies within the experimental error [TQ| ITT] . In 
order to clarify this, we have in this work studied the oxy- 
gen vacancy formation in Nao.75Co02 by first principles 
density functional theory (DFT) calculations, at various 
temperatures and oxygen partial pressures. We found 
that oxygen vacancies have high free energies of forma- 
tion even at high T, suggesting that the experimentally 
observed weight loss may have a different origin. 
We assumed the following reaction (using Kroger- Vink 
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Figure 1: Structure of the three Nao.rsCoCb models, in their 
starting (non relaxed) positions, as looking along the 001 di- 
rection. Top panel: Diamond, Bottom left panel: FHC, Bot- 
tom right panel: ZigZag model. Nal sites are shown as light 
blue balls, and Na2 as yellow balls. Co ions are shown as 
dark balls at the center of the blue Co06 octahedra. All the 
oxygen sites are shown in red, but the ones removed in the 
vacancy calculations, shown in green. 

notation): 

Na x Co0 2 v± Na x Co0 2 -5 + SVg + qSe' + -0 2 (1) 

The Gibb's free energies of formation for oxygen vacan- 
cies AG f(P,T) at a pressure P and temperature T were 
accordingly computed from 

AG f (P,T) = G def (P,T)-G perf (P,T)+^ 02 (P,T)+q» e 

(2) 

where q is the charge state per defect, and is the 
chemical potential of the species i. Thermal contribu- 
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tions of the solid were considered as negligible with re- 
spect to those of the oxygen gas, whose chemical potential 
jao 2 (Pj T) was calculated from 

^o 2 (P,T) = Ho 2 (P ,0) + ^o 2 (P°,T)+k B Tln^^j , 

(3) 

where ks is the Boltzmann constant. The dependence 
of no 2 on the temperature A/io 2 was taken from 
thermodynamic tables, while enthalpies at T = K and 
standard pressure were taken as the calculated DFT 
total electronic energies. System charge q was simulated 
by adjusting the total number of electrons in the cell 
and at the same time adding a compensating jellium 
background to avoid diverging Coulomb contributions. 
The electron chemical potential was determined by 
the Fermi energy corrected to account for the charge 
adjustment [12J. 

Periodic DFT as implemented in the VASP package [13J 
was used throughout this work. A spin polarized gradi- 
ent corrected Perdew-Burke-Ernzerhof (PBE) functional 
[T4] was used with a plane wave energy cutoff of 600 
eV. Core electrons were included through the projector 
augmented wave (PAW) method [15J. The reciprocal 
space was sampled by a T-centered k-point grid in which 
the maximum distance between points is 0.15 x27r/|a|. 
Ionic positions were relaxed until the maximum force 
was lower than 0.03 eVA -1 . Simultaneous relaxation of 
the lattice parameters was performed for the defect free 
structures. 

To sample the many possible sodium ordered structures 
we adopted three different models for Nao.75Co02. They 
are shown in Fig. [T] and named diamond, filled honey- 
comb (FHC) and zigzag. They are based on a tetragonal 
y/3 x 4 x 1 supercell consisting of 16 formula units [16J 
and they represent different occupation ratios (1, 1/2, 
and 1/5) of the two Na crystallographic inequivalent 
sites: Nal and Na2. The three models are based on 
the experimental structure reported by Zandbergen et 
al.[2\ (Diamond), and on the theoretical study of Meng 
et a/.[3l [17] (FHC and ZigZag). These were chosen 
as a representative selection of simple periodic models 
from the large manifold of structures, suitable for a first 
principles study. Nevertheless, Na self-diffusion is an 
efficient process already at room temperature [IB], so we 
expect a strong Na disorder at high T. 
We considered oxygen vacancies at three different 
positions for each of the three models, shown as green 
balls in Fig. [T] We first note that the models with 
oxygen vacancies are half- metallic (not shown here), 
similarly to the stoichiometric Nao.75Co02 structure. 
Oxygen vacancies were considered in their neutral, +1 
and +2 oxidation states, hence as Vq, V* and V**, 
and their free energies of formation (AGf) are shown 
in Table l[] for the high temperature regime (T = 1000 
K) relevant for applications. Results are reported for 
the three different models in both a y/3 x 4 x 1 and 
y/3 x 8 x 1 supercell, corresponding to Nao.75Co02-5 
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Figure 2: Formation energy for the doubly charged oxygen va- 
cancy as a function of the Na content. Each formation energy 
has been calculated from a few structures with different or- 
dering, mean values and standard deviations are shown. Line 
is a guidance for the eye. See text for details. 



AG f (eV) 
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vs # v v* 


V - V V* 


Diaml 


2.19 2.26 2.30 


2.27 2.27 2.28 


Diam2 


2.30 2.32 2.35 


2.32 2.34 2.36 


Diam3 


2.29 2.33 2.36 


2.35 2.37 2.39 


ZigZagl 


2.35 2.37 2.37 


2.36 2.37 2.39 


ZigZag2 


2.32 2.37 2.39 


2.39 2.40 2.40 


ZigZag3 


2.37 2.42 2.44 


2.41 2.43 2.45 


FHC1 


2.46 2.44 2.48 


2.40 2.39 2.42 


FHC2 


2.37 2.39 2.41 


2.38 2.37 2.39 


FHC3 


2.19 2.24 2.26 


2.20 2.23 2.25 



Table I: Defect formation free energy as defined in Eq. [2] 
for neutral, singly and doubly charged oxygen vacancy in 
Nao.75Co02-<5 computed for two different 5 values. T = 1000 
K and P(Q2) = 1.0 atm. All the values are in eV. 



with S = 0.031 and 5 = 0.016, respectively. The table 
values were based on the Fermi energy Ep of the perfect 
material and P(02) = 1.0 atm. 

Formation energies are distributed in a narrow range of 
0.2 eV, with no clear preference for any of the structural 
models. This suggest that the Na ordering only plays 
a minor role for the vacancy formation energy. The 
distribution of formation energies is also narrow for 
the two vacancy concentrations 5; lower than 0.1 eV 
in most of the cases. This means that the Coulomb 
potential generated by the charged defect is efficiently 
screened already at the length scale of the shortest 
lattice parameter, i.e. ~5.0 A. Nevertheless, it appears 
that the doubly charged defect V** is slightly favored 
compared to the other two possible charges, suggesting 
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Figure 3: Formation free energy (left) and concentration 
(right) of doubly charged oxygen vacancies as a function of 
oxygen partial pressure. Three different temperatures are in- 
vestigated. 



a non-perfect dielectric screening (Fermi liquid). 
Overall, the defect formation energies shown in Table 
l[] are quite high, ranging from 2.19 to 2.48 eV at 
T = 1000 K and P(02) = 1.0 atm, implying a very low 
concentration of vacancies. Since the Na ions are very 
mobile at high temperature, one could expect large local 
variations of the sodium concentration, modifying the 
local vacancy formation energies. We therefore studied 
the Vq" defect formation for a range of Na concentrations 
x. Two structures were generated for each x value, and 
the oxygen vacancy was created at three different sites 
for each Na arrangement. Averaged results are shown in 
Fig|2j using the same conditions as in Table l[] When the 
Na content is reduced, AG/ lowers almost linearly down 
to £=0.50, where it reaches a plateau of about 1.70 eV. 
The change in the formation energies is substantial, but 
in Nao.75Co02 it is unlikely that diffusion alone could 
create (even locally) low enough Na concentrations to 
reduce AG/ by e.g. 1 eV. 

In order to quantify the vacancies concentration under 
experimental conditions we calculated AG/ for x = 0.75 
at different temperatures for the structure with lowest 



free energy of formation (2.18 eV) in Table l[] The 
results are shown in the left panel of Figure [3] for three 
different temperatures (700, 1000 and 1300 K) as a 
function of P(C>2). In the right panel we show the 
corresponding vacancy concentration per formula unit, 
S. We also note that our formation energies are in good 
agreement with those recently reported by Yoshiya et 
al. [19J for x = 0.50 and x = 1. It is clear that the 
oxygen vacancy free energy of formation is very high 
at relevant temperatures and oxygen partial pressures. 
This conclusion holds also when taking into account 
errors arising from the methodology or approximations 
we adopted. Even at extreme experimental conditions 
such as T = 1300 K (close to the homogeneous melting 
point[20j) and P(0 2 ) = 10~ 7 atm the vacancy con- 
centration would be as little as S ~ 6.5 x 10 -4 . If 
more common experimental conditions are chosen, e.g. 
T = 700 K and P(0 2 ) = 10" 1 atm, we obtain S ~ lO" 18 . 
This is several orders of magnitude lower than reported 
in some previous publications, e.g. Ref. |6J. 
Our results may help to settle the disagreement between 
experimental reports on the oxygen stoichiometry in 
Na ;E Co02. At high sodium content, concentrations of 
vacancies ranging from 5 = 0.05 to 0.16 have been 
reported [6j [8[ |9], but other studies have concluded that 
the oxygen concentration is indeed stoichiometric [T0| [TT] . 
Our results are clearly consistent with the latter ones. 
We can here only speculate on the reason for this 
discrepancy. One possible explanation is existence of 
secondary phases in some of the experiments, where 
these phases are responsible for the oxygen weight 
loss upon reduction of oxygen partial pressure. Other 
possibilities comprise Na evaporation, the presence 
of other defects than considered here, or unintended 
interaction between the sample and experimental setup. 
In conclusion, we showed that the equilibrium concentra- 
tion of oxygen vacancies in sodium cobaltate is very low 
at conditions relevant for applications as thermoelectric 
material, solving a long lasting experimental controversy. 
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